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Recent  calculations  [S.  Prager,  AFOSR-2803  (June  1962)]  of 
upper  bounds  for  the  effective  diffusion  coefficient  (or  conducti¬ 
vity)  in  porous  media,  in  terms  of  certain  statistical  parameters 
of  the  random  geometry,  are  reformulated  so  as  to  apply  speci¬ 
fically  to  a  bed  of  spherical  particles.  The  calculations  are 
simplified  by  considering  an  idealized  bed  in  which  centers  are 
randomly  situated  without  restricting  the  spheres  to  nonoverlap¬ 
ping  locations.  The  result,  applicable  to  randomly  overlapping 
spheres  of  either  uniform  or  nonuniform  sizes,  gives  the  upper 

bound  for  the  effective  diffusion  coefficient  as  4>D  /[l  -  —  fn<)>], 

2 

where  Dq  is  the  actual  diffusion  coefficient  in  a  fluid  which  fills 
the  void  regions  of  the  bed  and  ]>  is  the  void  fraction.  This 
result  is  compared  with  experimental  results  by  various  investi¬ 
gators  for  nonoverlapping  spheres  and  also  with  the  best  upper 
bound  that  can  be  calculated  without  taking  the  statistics  of  a 
particular  random  geometry  into  account  [Z.  Ha  shin 
S.  Shtrikman,  J.  Appl.  Phys.  33,  3125  (1962)]. 
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1.  INTRODUCTION 

Attempts  at  estimating  the  extent  to  which  the  apparent  diffusion 
coefficient  is  lowered  when  diffusion  occurs  in  the  open  regions  of  a 
porous  medium,  and  analogous  calculations  related  to  the  thermal  or 
electrical  conductivity  of  inhomogeneous  substances,  are  at  least  as 
old  as  the  writings  of  Maxwell,  *  The  purpose  of  the  present  work  is 
to  show  how  previous  theoretical  estimates  can  be  improved  by  the 
utilization  of  statistical  information  about  the  random  geometry  of  the 
medium. 

Maxwell's  formula*  gives  the  electrical  conductivity  of  a  medium 
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which  contains  a  dilute  suspension  of  spheres;  as  applied  to  nonconduct¬ 
ing  spheres  the  result,  expressed  in  the  terminology  of  our  diffusion 
problem,  is 

De/Do  =  4>/[l  +i(l  -<►)].  (1) 

where  4>  is  the  volume  fraction  not  occupied  by  spheres,  Dq  is  the 
diffusion  coefficient  in  the  open  space,  and  is  the  effective  diffusion 
coefficient  defined  as  the  ratio  of  average  flux  to  average  concentration 
gradient.  It  has  recently  been  shown  by  Hashin  and  Shtrikman^  that 

the  right  side  of  this  equation,  originally  derived  for  <$> - ^  1,  is  an 

upper  bound  for  D^/Dq  in  an  isotropic  medium  for  any  value  of  <J>  even 
when  the  solid  portions  of  the  medium  are  not  spheres.  Moreover 
these  authors  showed  that  this  bound  is  the  closest  that  can  be  obtained 
without  providing  additional  statistical  specification  of  the  medium, 
t  A  method  for  introducing  the  required  statistics  has  recently 
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been  presented  by  Prager  in  terms  of  a  variational  formulation  of  the 
diffusion  problem.  His  results  are  applicable  to  a  very  general  class 
of  isotropic  porous  media;  however,  the  statistical  data  required  for 


3. 


such  applications  are  not  yet  available.  Therefore,  in  the  present  work, 
the  variational  approach  is  reformulated  so  as  to  apply  specifically  to  a 
bed  of  spheres.  It  will  be  seen  that  the  appropriate  statistical  specifica¬ 
tion  can  then  be  greatly  simplified. 


2.  EFFECTIVE  DIFFUSION  COEFFICIENT 

The  bound  on  O  I D  ,  which  is  formulated  in  reference  3  in  terms 
e  o 

of  a  trial  concentration  fluctuation  c',  is  equivalent  to  the  following 
inequality  in  terms  of  a  trial  concentration  c: 
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where  g(r)  is  a  random  function  of  the  position  vector  r,  defined  to 

A/y  W 

have  the  value  zero  in  the  solid  regions  and  unity  in  the  void  regions  of 
the  porous  medium.  The  angular  brackets  denote  an  average  taken  over 
a  volume  V  which  is  large  in  comparison  to  the  scale  of  the  inhomo- 
geneties  of  the  medium.  For  example,  the  "porosity"  or  void  fraction 
(j>  is  given  by 


(3) 


2.  1  Trial  Function  for  a  Bed  of  Spheres 
The  trial  concentration  gradient  which  is  now  used  in  Eq.  (2) 
involves  a  sum  of  contributions  from  each  sphere  center,  the  i**1  center 
being  located  at  the  position  r..  To  this  sum  is  added  a  constant 
vector  A  in  the  direction  of  the  average  concentration  gradient  <^c^.. 
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r  -  r.  and,  for  spheres  of  radius  a 
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This  form  of  the  contribution  from  each  center  is  suggested  by  the 
solution  of  the  diffusion  problem  in  which  the  obstruction  is  a  single 
isolated  sphere}  *  its  use  here  insures  that  the  upper  bound  for  D 

c 

becomes  exact  as  4>  approaches  unity. 

The  scalar  multiplier  m  in  Eq.  (4a)  is  to  be  adjusted  so  as  to 
minimize  the  right  side  of  Eq.  (2)  when  substitution  is  made  for  <^c 
from  (4).  The  result  of  this  substitution,  after  some  reduction  involving 


use  of  Eq.  (3),  can  be  written  as 
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fe  <  »A  +m  \g(#t\Vf(4i'i)]> 
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where  n  is  the  average  density  of  sphere  centers. 

The  optimum  value  of  m,  obtained  by  differentiating  the  right 
side  of  Eq.  (5)  with  respect  to  m,  is 
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Expanding  the  squared  sum  in  the  denominator  and  substituting  the 
resulting  expression  for  m  into  Eq.  (5)  gives 

<t>  /  q  (6a) 
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The  parameter  q  is  related  to  a  "tortuosity  factor"  which  has  been 


(6b) 


discussed  by  Carman  and  others. 
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2,2  S  irnplif  .i.CAt  J  o.r>  of  the  Random  G  a  ora  e  try 

It  is  in  the  evaluation  of  the  summedions  over  sphere  centers  and 

over  pairs  of  centers,  which  appear  in  Eq.  (6b),  that  the  statistical 

properties  of  the  bed  of  spheres  must  be  considered.  In  this  connection, 

one  naturally  thinks  of  a  bed  in  which  each  center  is  at  least  one  diameter 

from  its  neighbors,  but  the  evaluation  of  the  pertinent  sums  is  then  far 

from  simple.  However,  it  is  possible,  at  least  conceptually,  to  produce 

a  bed  of  spheres  for  which  the  calculations  are  greatly  simplified;  this 
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is  achieved  by  removing  the  constraint  on  the  locations  of  centers. 

Thus,  if  centers  are  distributed  at  random,  regardless  of  inter¬ 
penetration  of  spheres,  then  the  coordinates  p  .  and  p.  of  the  members 

/  va  1  %ia  ^ 

of  an  ij  pair  of  centers  in  Eq.  (6b)  are  completely  uncorrelated  and  it  is 

seen  from  the  form,  Eq.  (4b),  of  Up.)  and  f(  p.)  that  the  double  sum  in 

/  1  1  J 

Ws  ** 

(6b)  vanishes. 

The  single  sum  in  Eq.  (6b)  is  also  readily  evaluated  when  the 
spheres  are  allowed  to  overlap;  for  then  the  probability  that  a  random 
point  is  in  a  void  region,  even  though  the  point  is  known  to  be  in  the 
vicinity  of  the  i  center,  is  simply  $  for  all^CK  >  a.  Thus,  when 

IVY\ 

Eq.  (4b)  is  used  for  f(p  the  single  summation  over,  say,  N  spheres 

Art 

can  be  expressed  as 
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To  eliminate  n  and  a  from  the  final  result,  we  note  that 

4>  =  exp  (-4ira^n/3);  (8) 

this  expression  for  the  void  fraction  of  a  bed  of  randomly  overlapping 
spheres  is  derived  in  the  Appendix  along  with  certain  other  interesting 
statistical  properties  of  this  kind  of  porous  medium. 
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Substitution  of  Eq.  (7)  in  (6b)  and  then  using  (8)  gives 

q1  =  1  -  (9a) 

2.3  Nonuniform  Spheres 

Although  Eq.  (9a)  was  derived  for  a  bed  of  uniform  spheres,  it 
is  interesting  to  note  that  the  same  result  is  obtained  when  the  bed 
contains  overlapping  spheres  of  several  sizes.  The  generalized 
derivation  is  a  straightforward  repetition  of  the  procedure  used  above 
for  uniform  spheres;  we  need  only  mention  here  the  form  of  the  trial 
function  and  the  expression  for  the  void  fraction. 

The  trial  function  which  replaces  Eq.  (4)  is 


V  c(r)  =  A  +  mj  \7fx(  _•.)  +  m2  +  —  (10a) 
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Thus  there  is  an  adjustable  multiplier  im  for  each  different  sphere 

radius  a.. 

) 

The  required  expression  for  the  void  fraction,  calculated  by 
the  same  method  as  is  described  for  uniform  spheres  in  the  Appendix,  is 


$  =  exp  [-(4tt/3)  (n^j3  +  n2a23  +  ---)  ],  (11) 

where  n^,  n2, - are  the  respective  densities  of  sphere  centers  for 

the  several  sizes. 


3.  DISCUSSION  OF  RESULTS 

The  result  given  by  Eq.  (9a)  may  be  compared  with  the  corres¬ 
ponding  expression  for  q  from  Maxwell's  formula,  Eq.  (1): 

v  -  1  +  \  <> 


(9b) 
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We  note,  for  example,  that  for  dilute  beds  (4>  — ^  1)  the  two  expres¬ 

sions  are  identical  to  first  order  in  (1  -  4>).  A  graphical  comparison 
of  (9a)  and  (9b)  is  shown  in  Fig.  1. 

The  values  q.j  and  qm  may  be  regarded  as  estimates  q  of  the 
true  value  q*  defined  such  that  Dg/Do  s  <|>/q*.  Then,  from  (6a),  these 
estimates  are  rigorous  lower  bounds  for  q*  and  the  previously  cited 
theorem  of  Hashin  and  Shtrikman^  asserts  that  qm  is  the  greatest 
lower  bound  that  can  be  derived  without  considering  statistical  descrip¬ 
tions  of  the  medium  other  than  the  void  fraction.  That  a  higher  bound 
is  indeed  obtained  when  additional  statistics  are  introduced  is  illustrated 
by  the  graph  of  q^  vs.  <f>,  shown  as  the  upper  curve  in  Fig.  1. 

It  is,  of  course,  expected  that  use  of  sufficient  statistical  infor¬ 
mation  would  provide  bounds  on  q*  which  are  sufficiently  close  to 
furnish  useful  predictions  of  experimental  results.  Although  the  kind  of 
sphere  bed  for  which  Eq.  (9a)  provides  a  rigorous  bound  would  be 
difficult  if  not  impossible  to  produce  in  the  laboratory,  it  is  nevertheless 
interesting  to  compare  calculated  values  of  q^  with  measured  values  of 
q*  as  obtained  from  more  readily  fabricated  sphere  beds  in  which 

centers  are  separated  by  at  least  one  diameter.  Such  measured  values 
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by  various  investigators  have  been  tabulated  by  Carman.  These  results 
were  obtained  for  various  kinds  of  porous  media  among  which  were  beds 
of  uniform  spheres  and  also  mixtures  of  two  and  three  different  sphere 
sizes;  the  data  for  the  nonuniform  as  well  as  uniform  sphere  beds  are 
plotted  in  Fig.  1.  This  comparison  of  experimental  data  for  nonover¬ 
lapping  spheres  with  the  theoretical  estimate  for  overlapping  spheres 
(upper  curve)  cannot,  of  course,  be  regarded  as  a  valid  test  of  the  latter; 
it  does  provide,  however,  a  good  indication  that  the  kind  of  rigorous 
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bounds  calculated  here  can  indeed  furnish  useful  estimates  of  effective 
diffusion  coefficients  or  conductivities. 
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APPENDIX:  VOID  FRACTION  AND  TWO-  AND  THREE - 
POINT  CORRELATIONS  FOR  A  BED  OF  OVERLAPPING  SPHERES. 


To  obtain  an  expression  for  the  void  fraction  <|>  in  terms  of  the 
density  of  sphere  centers  n  and  the  sphere  radius  a,  we  first  calculate 
the  probability  of  placing  N  centers  at  random  in  a  large  finite  volume  V 
in  such  a  way  that  a  smaller  volume  v  contains  no  centers.  Since  each 
random  placement  of  a  center  is  independent  of  the  positions  of  the  other 
centers,  this  probability  is  [(V-v)/V]^  =  [1  -  (nv/N)]^.  As  N  is  made 
larger  and  larger,  holding  n  and  v  fixed,  the  probability  Py  that  the 


volume  v  contains  no  centers  becomes 


lim 


P  = 
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N-> 


(1  - 


nv 


N 


N 

)  =  e 


-nv 


(Al) 


Now  4>  can  be  regarded  as  the  probability  that  a  random  point  is 
not  contained  by  any  sphere  in  the  bed,  or,  equivalently,  as  the  proba¬ 
bility  that  the  portion  of  the  bed  which  lies  within  a  distance  a  from  the 

random  point  contains  no  centers.  Hence  <j>  is  given  by  Eq.  (Al)  with 
3 

v  =  4ira  / 3 ,  which  is  the  result  given  in  Eq.  (8). 
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It  is  interesting  to  note  that  Eq.  (Al)  also  provides  an  expression 

for  the  two-point  correlation  function  S(s^)  defined  as  the  probability  that 

two  points  separated  by  a  distance  r  both  lie  in  void  regions.  This 

or  closely  related  functions  have  been  used  in  the  previously-cited 

3 

calculation  of  diffusion  in  porous  media  and  in  calculations  of  x-ray 

scattering  in  random  media^  and  of  viscous  flow  of  fluids  through  porous 
?  8 

media.  *  It  can  be  obtained  for  a  bed  of  overlapping  spheres  by  setting 
v  in  Eq.  (Al)  equal  to  V,  .  (r)  as  defined  in  reference  8.  The  result  is 


S(r)  = 


-  r  4-rra  n  ..Sr 
exp  [-  - - (1  +  — 

3  4a 


(^exp  [• 
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8ira  n 


r  <  2a 


r  >  2a 


(A2) 


This  form  of  the  two-point  correlation,  since  it  is  derived  for  a  specified 
(albeit  artificial)  random  geometry,  may  be  preferable  to  the  simpler 
exponential  form  discussed  by  Debye,  Anderson,  and  Brumberger;  ^  it 
is  not  certain  whether  there  is  any  three-dimensional  geometry  which  is 
consistent  with  the  kind  of  randomness  postulated  by  these  authors. 

The  three -point  correlation  function  G(r,  r')  has  also  been  used 

•*H  /Vii 
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in  previous  calculations  for  porous  media.  ’  ’  It  is  obtained  by  setting 

v  in  (Al)  equal  to  Vv  ,  (r,  r')  as  defined  in  reference  8.  The  latter  quantity 

Mr-,  m\ 

represents  the  volume  enclosed  by  three  spheres  with  centers  at  the 

vertices  of  a  triangle  with  sides  r,  r',  and  r  -  r';  the  cumbersome 

-V,  s-A  W  ,v^ 

expression  for  this  volume,  involving  various  forms  for  various  kinds 
of  r,  r'  configurations,  will  not  be  repeated  here. 
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FIGURE  CAPTION 


Fig.  1,  Tortuosity  parameter  q  =  <j>D^D  vs.  void  fraction  <j> 
for  a  bed  of  overlapping  spheres.  Upper  curve  -  eq.  (9a),  lower 
curve  -  eq.  (9b).  The  plotted  points  represent  experimental  results 
for  nonoverlapping  spheres  from  reference  4:  ®  -  bed  of  uniform 
spheres;  O  -  mixtures  of  spheres  of  more  than  one  size. 

Both  curves  are  theoretical  lower  bounds;  the  upper  curve 
illustrates  the  improvement  obtained  by  incorporating  in  the  calcula¬ 
tion  a  statistical  description  of  the  porous  medium. 
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